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Abstract 

Overlap fermions have an exact chiral symmetry on the lattice and are thus an appropriate tool 
for investigating the chiral and topological structure of the QCD vacuum. We study various chiral 
and topological aspects of quenched gauge field configurations. This includes the localization and 
chiral properties of the eigenmodes, the local structure of the ultraviolet filtered field strength 
tensor, as well as the structure of topological charge fluctuations. We conclude that the vacuum 
has a multifractal structure. 
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I. INTRODUCTION 

Overlap fermions possess an exact chiral symmetry on the lattice [3] and realize the 
Atiyah-Singer index theorem at finite cut-off 
local definition of the topological charge density 
a powerful tool for investigating the chiral and topological structure of the QCD vacuum. 



. Furthermore, they lend themselves to a 

n 

5(. Altogether, this makes overlap fermions 



In this paper we shall address the vacuum o 



Preliminary results have been published in 



the quenched theory at zero temperature. 



B, q y, [9 



9]. In a forthcoming paper we shall 



extend the investigation to dynamical QCD at zero and finite temperature. First results 
obtained for full QCD in the vicinity of the chiral phase transition have been reported in 

IK 



It is known for some time 
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12j that the long distance properties of QCD are well 
described by the low-lying eigenmodes of the overlap operator. The question of low-mode 
dominance has been raised earlier in Ref. [13I . [l^ . In Fig. [T] we compare the full pion 
propagator with the truncated one spanned only by the 40 lowest eigenmodes. 

We see that both propagators tend to each other at distances > 0.5 fm. A suitable 
truncation of the overlap operator thus acts as an infrared filter, which allows us to sepa- 
rate the truly nonperturbative degrees of freedom of the QCD vacuum from the ultraviolet 
noise. This is in accord with conventional wisdom, namely that chiral symmetry breaking is 
encoded in the low-lying modes of the Dirac operator and both, chiral symmetry breaking 
and confinement have the same dynamical origin. To make this more precise, the vacuum 
structure, seen changing with an increasing number of fermionic modes, is the general theme 
of this paper. Thus, we shall come back also to the early saturation of the pion propagator 
in more detail later in this paper. 



Earlier lattice investigations of QCD vacuum structure as reviewed in Refs 
partly relied on gauge fixing and subsequent projection onto appropriate subgroups of color 
SU(3). The topics under discussion there are mainly related to the confinement of heavy 
quarks. In contrast, with the gauge invariant overlap approach, we come closer to the 
confinement issue of light quarks. As a particular advantage, it establishes a direct link 
between topological excitations and light quark propagation, which was missing so far. After 
all, we believe that light fermions are a major element of the low-energy effective action of 
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FIG. 1: The full pion propagator compared with the truncated pion propagator constructed as 
an all-to-all correlator summing over the 40 lowest modes only with a quark mass Tp = 0.01 (see 
Eq. (fT2l) ). This calculation uses 250 configurations on the 16 3 x 32 lattice generated at j3 = 8.45. 



QCD 



id- 



In the instanton model (as a prototype of other semiclassically motivated models), with 
the almost-zero-mode band approximation, this link is realized in a very economic way. We 
summarize these models by saying that they all are based on classical, selfdual or antiself- 
dual solutions of the Euclidean equations of motion carrying zero modes of definite chirality. 
While the vacuum fields are formed in these models as superpositions of such extended 
solutions, in the almost-zero-mode band approximation the spectrum of the lowest quark 
eigenstates in the complex vacuum is analogously formed by linear combinations of the cor- 



responding zero modes. Examples are the instanton liquid model (for a review see Ref. 19j ) 
where this program has been carried through, or the caloron gas model [2o| where this has 
still to be done. Both models are based on gauge fields piecewise coherent over O(0.5 fm) 
but with decaying field strength correlations beyond that distance. In order to justify these 
models, a good deal of lattice vacuum studies has been devoted to instanton and caloron 
searches. 



In the past, cooling [2l|, [22], underrelaxed cooling [231] , restricted improved cooling [24] 
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and other smoothing techniques including inverse blocking [25|, |26|, |27|, renormalization 



group cycling 



281 ]. finally rationalized as smearing [29|], have been used for this purpose. 



Since these methods contain a bias towards classical solutions, the results have been taken 
with certain reservations. In particular, the number density of lumps ("instantons and 
antiinstantons" ) was strongly influenced by systematical effects [30[]. Only recently, when 
valence overlap fermions have become a viable option, a local definition of the topological 
charge density, q(x) oc Tr (F llv F jt J), in terms of the overlap Dirac operator (see Eq. (fT8l) ) 
became practicable (see e.g. Ref. [31]). Surprisingly, for Monte Carlo configurations this 
density shows a three-dimensional, laminar and highly singular structure which seems to 
rule out the instanton and similar pictures. 

Easier accessible are mode-truncated, i.e. ultraviolet filtered versions of the topological 
charge density (see Eq. ([15]) ). Still with large computational effort, such evaluations of 
the topological charge density are now superseding cooling or smearing techniques. The 
latter techniques were necessary to be applied in the past to the gauge field before the 
topological charge density could be calculated according to a gluonic definition (see for 
example Ref. 32]). The number of fermionic modes (ultraviolet filtering) now replaces the 
number of smearing iterations as the control parameter of smoothness. The overlap definition 
of the topological charge density, if recast [3l| into its spectral representation (Tl9l) . can be 
evaluated with any desired degree of ultraviolet smearing represented by a cut-off A cut . In this 
paper the two aspects (the infrared and the ultraviolet one) of the topological structure will 



be described. Whereas the infrared structure is definitely associated with chira 



breaking, the relation of the surprising divergent structure in the ultraviolet 
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symmetry 



33 



34| to 



351 ] . is still hypothetic. 



the phenomenology of confinement, as stressed by 

The all-scale topological charge density ( 1181) . containing fluctuations of all scales, from 
the lattice spacing a to global structures percolating through the full lattice volume, is 



calculable 



311 ] directly in terms of the overlap operator. This does not require to know its full 



spectrum and all eigenmodes. It is calculable, although computationally very demanding, 
as a local trace due to the form of the Neuberger overlap Dirac operator. In contrast to this, 
the ultraviolet filtered topological charge density can be quite easily calculated. Moreover, 
the eigenmodes can be used for many applications , ISj. In this paper, for instance, 

we will elucidate the low-mode dominance of the pion propagator in more detail. 
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In the past, selfdual objects have 
coherent lumps of topological charge 



Deen searchec 
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'or on the 
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26|, 
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attice primarily by looking for 
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29j which play a prominent 



role in models of vacuum structure. When the overlap-based topological charge became 



36. 



371 ] . Early arguments [35J against 



calculable, these models were heavily challenged 
the instanton dominance of the functional integral have been raised again. The striking new 
argument against the instanton structure was that the four- dimensional extendedness typical 
for semiclassical lumps cannot be reconciled with the newly discovered three-dimensional 
sign-coherent global structures 33j which are infinitely thin in the co-dimension (equal to 
one 34J]). This is the picture that has emerged for the all-scale topological density. All what 
can be said in defense of instantons etc. is that it is unknown how the effect of quantum 
fluctuations modifies the topological charge density profile of a classical instanton or caloron 
which is described by overlap fermions in full agreement with the classical profile. 

In this paper, concerned with fully quantum lattice configurations, we will further develop 
this picture. It is pointed out that the topological charge density, in regions of higher density, 
possesses lower- dimensional connected structures. In the full four-dimensional landscape we 
find all substructures, from zero-dimensional peaks over one-dimensional ridges, to two- 
dimensional walls. They are all nested inside the three-dimensional sign coherent structures 
mentioned above. On the other hand, it will become clear in which sense a model based on 
a dilute gas of (anti-) selfdual domains can also be supported by an overlap- fermion based 
vidently, for this an ultraviolet filtered version of the overlap-based field strength 
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401 ] is necessary, supplemented by the mode-truncated topological charge 



311 ] . Both can be projected out by a certain small number of true low-lying overlap 



analysis, 
tensor 
density 
modes. 

This latter point of view is qualitatively attractive because the low-lying fermionic modes 
are known to encode the most important phenomenological properties of the QCD vacuum 
relevant for the physics of light hadrons. The case for concentrating on the lowest part 
of the Dirac spectrum is illustrated by the fact already mentioned that a modest num- 
ber of eigenmodes is sufficient to reproduce the propagator of the lightest (pseudoscalar) 
[13I . Q (see also Fig. [T]). The lowest-mode dominance allows for a system- 



hadrons 
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atic improvement of the propagator by the so-called low-mode averaging method 41 
In this and a following paper this will be investigated more carefully. 
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43]. 
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Recently, the localization properties 



451,1461,147, 



36 



48 Land t 



re effective dimension of 



23, 
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50J have received strong 



the lower eigenmodes, and their local chiral properties 
interest. The results partly contradict each other. Concerning the localization the situ- 
ation has been reviewed in Ref. 5l|. There, de Forcrand points out the basic difference 
between SU(2) fields simulated with Wilson action, when scale a dislocations have too low 



action 



52] to compensate for their entropy ~ log a 4 , and SU (3) gauge fields simulated with 



Luscher-Weisz action with dislocations suppressed. The general motivation behind this high 
attention for issues like localization and (fractal) dimensionality was the hope 53|, |54| that 



the lowest modes are bound to certain singularities of the gauge field that in turn would 
explain confinement. Here, "singular" is understood in a sense opposite to semi classical 



lumps, just in the spirit of Witten's 
more recent point of view in Refs 
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criticism of instantons. This is discussed from a 
56] . Uncorrelated instantons, indeed, are found unre- 
lated to confinement. It should be mentioned in passing, however, that other semi classical 
configurations with non-trivial asy mptotic holonomy can be related to confinement, as was 
recently demonstrated in Ref. 20]. 

Another source of inspiration was the hypothetical analogy drawn between the finite^ 
temperature chiral transition to the Anderson transition (metal insulator transition) 53, 



58, 



59] in condensed matter physics. As precursor of the chiral symmetry restoring transition, 



before the spectral gap opens, the spectrum in the gap region is expected to become critical, 
exposing both spectral criticality and multifractality of modes. We will come back to these 
aspects in a forthcoming publication dealing with finite temperature full QCD. In the present 
case of quenched QCD the whole low-lying spectrum seems to be critical and multifractal. 

This paper is organized as follows. In Section [Til the basic tools and settings of the 
quenched simulation and for handling the eigenmodes of the overlap operator are explained. 
In Section II III the properties of individual modes are discussed. This begins with the distri- 
bution of the topological charges Q (given by the number and chirality of zero modes) and 
the spectral density of non-zero modes. From the latter, the chiral condensate is extracted 
by fitting our anrte-voaune data to jetted data, perturbation theory Q. Then the 
localization behavior, the dimensionality and the local chiral properties (chirality X(x)) of 
individual modes are discussed. In Section [TV] we describe our experience with the topo- 
logical charge density defined Q through the overlap Dirac operator and its ultraviolet 
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filtered (mode-truncated) variant. The existence of a three-dimensional, singular and sign- 
coherent global structure revealed by the all-scale topological charge density, first pointed 
out in Ref. 33], is confirmed and described in more detail. The negativity of the two-point 

n n 

function of the topological charge density [6J, 161(1 is found to be realized for the all-scale 
topological charge density for high enough j3, i.e. at lattice spacing a < 0.1 fm. In the same 
Section HVl we also discuss the topological structure in terms of what we call g-clusters. The 
multiplicity, size, distance and the percolation behavior of g-clusters is compared for various 
levels of ultraviolet filtering (A cut ) and for the all-scale topological charge density. With an 
appropriate A cut the g-cluster structure of the mode-truncated density gA cut ( x ) can be made 
agreeing with the .R-cluster structure defined later in Section |Vj Finally, estimating the 
fractal dimension of g-clusters related to the all-scale topological charge density at various 
threshold values g cut , we find that the three-dimensional structure of sign-coherent topo- 
logical charge is supplemented by nested, lower- dimensional substructures at higher density 
|g(x)|. Obviously these features have no analog for the mode-truncated topological charge 
density. In Section |V] an ultraviolet filtering technique for the field strength tensor is in- 
troduced. It was inspired by a similar work by Gattringer 3^|. The representation of the 
field strength tensor in terms of the overlap Dirac operator was later discussed by Liu et 
al. 0, . An infrared field strength can be obtained that allows one to assign at each point 
in space-time the degree R(x) of selfduality or antiselfduality. Connected clusters (called 
i?-clusters) are found such that in the interior (anti-)selfduality is approximately satisfied. 
These clusters start to percolate and form a dilute network if slight deviations from (anti- 
selfduality are tolerated. Thus, locally perfectly self dual and antiselfdual vacuum fields 
are embedded in a vacuum that is globally neither selfdual nor antiselfdual. In Section IVII 
we summarize our findings and discuss from this point of view the saturation of the pion 
propagator. We draw conclusions and point out routes for further research. 

II. BASIC TOOLS 

A. Lattice ensembles 



Topological studies using the Wilson one-plaquette gauge field action suffer from dislo- 
cations [62J and should be treated with caution. For this paper the Liischer-Weisz gauge 
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field action [63] is used which is known to suppress dislocations and to greatly reduce the 
number of unphysical zero modes of the Wilson-Dirac operator. The Luscher-Weisz gauge 
field action is given by 



S[U]=f3 ~ ^ Tr (1 - f/ plaquette ) 

o 

plaquette 

+PrY1 ^ReTr(l-f/ rcctanglc ) (1) 

rectangle 

+PpJ2 ^ReTr(l-C/ p 



— v _ J parallelogram ) 

parallelogram 

in terms of the standard plaquette, the planar rectangle and the parallelogram loop terms. 
The latter are closed along the diagonally opposite parallel links on the surface of a 3D cube. 

The coefficients (3r and j3p are taken from tadpole improved perturbation theory 
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64 



65. 



P*. = ~Jh (1 + 0.4805 «) , (2) 

(3 P = -4 0.03325 a , (3) 
ul 

where 



no = ^- ( Re Tr /7 plaqucttc ) j , a = . (4) 

The couplings (5r and j3p have been selfconsistently determined through the calculation of 



the average plaquette in Ref. 



for a set of (3 values. We can write f3 = cop-, (3r = Cip? 
and (3p = c 2 p- with c + 8ci + 8c 2 = 1, which fixes the relation between (3 and the bare 

,2 



coupling g 

To investigate the volume dependence of our data we have simulated at three different 
volumes at fixed coupling (5 = 8.45. To explore the a dependence of the results also a 12 3 x 24 
lattice at (3 = 8.10 has been employed with approximately the same physical volume as the 
16 3 x 32 lattice at /3 — 8.45. Finally, a large ensemble of rather coarse configurations on a 
16 3 x 32 lattice at (3 — 8.00 became available in the QCDSF collaboration. A first physics 



analysis on this basis has been presented in Ref. [43j. The physical volume of this ensemble 
almost equals the biggest physical volume studied on 24 3 x 48 at (3 = 8.45. 

In Table I the statistics of lattices used in our investigation is listed. The lattice spacing 



for (3 = 8.00 and (3 = 8.45 was determined in 



43| , where we used the pion decay constant of 



R 
P 


\r 1 

CL [IUIJ 


T 3 v T 


o L s [fin] 


v [ini J 


Xtop 


~jf of config. 


jf of modes 


8.45 


0.105(2) 


12 3 x 24 


1.3 


5 


[167(3)MeV] 4 


437 


O(50) 


8.45 


0.105(2) 


16 3 x 32 


1.7 


16 


[169(3)MeV] 4 


400 


O(150) 


8.45 


0.105(2) 


24 3 x 48 


2.5 


81 


[168(4)MeV] 4 


250 


O(150) 


8.10 


0.142(2) 


12 3 x 24 


1.7 


15 


[171(l)MeV] 4 


251 


O(150) 


8.00 


0.157(3) 


16 3 x 32 


2.5 


74 


[172(4)MeV] 4 


2156 


O(170) 



TABLE I: Details of the quenched ensembles used in this study: couplings (3, lattice spacings a 
determined from the pion decay constant, lattice sizes (L s and Lt), physical lattice sizes (a L s 
and the volume V), topological susceptibilities Xtop> the statistics of configurations and available 
overlap fermion modes. 

f n = 92.4 MeV to set the scale. The value of the lattice spacing at (3 — 8.10 is interpolated 
jy fitting our measured pion decay constants at (3 = 8.00 and (3 = 8.45 to equation (6) of 
SSI- 
It is in this scale that the spatial linear extent a L s of the lattices, the 4D volumes V 
and the topological susceptibilities Xto P are obtained. The topological susceptibilities will 
be discussed in Section UlI Al 



B. Implementation of Neuberger overlap fermions 



Overlap fermions jl, 2] have an exact chiral symmetry on the lattice [3] and provide the 
cleanest known theoretical description of lattice fermions. Their implementation of chiral 
symmetry and the possibility to exactly define the index theorem on the lattice at finite 
lattice spacing allow one to investigate the relationship of topological properties of gauge 
fields and the dynamics of fermions. A further advantage of overlap fermions, in contrast to 
Wilson fermions, is that they are automatically 0(a) improved 69]. 
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The massless Neuberger overlap operator is defined by 



D(0) = P[l+'*L = ) t D W = M- P -, (5) 



a 



D W D W 



where we use the Wilson-Dirac operator Dw as input. M is Wilson's hopping term with 
r = 1. The negative mass parameter p is chosen to be 1.4, which represents a reasonable 





compromise between the physical requirement of good locality properties [43j, [70] of the 
overlap operator and a performance requirement, demanding a small condition number of 
H^y, where H\y = is the Hermitean Wilson-Dirac operator. In Fig. [2] the effective 

range of -D(O), represented by the decay of 

F(r) = /^max\D(0;x,y)\ ) ) , (6) 

with respect to the Euclidean distance 

M=(1^J , (7) 

is shown, as obtained on a 16 3 x 32 lattice at (3 = 8.45 for p = 1.4 , 

Asymptotically, F(r) oc exp (—p r/a), where p depends on p (and a priori also on the 
roughness of the configurations, i.e. 0). It turns out, however, that the fitted slope p = 
1.11(1) is practically independent of (3. Thus, the overlap Dirac operator is not ultralocal. 
It has a range that shrinks together with the lattice spacing towards the continuum limit. 
This is sufficient for any Dirac operator D to be called local. 

To compute the sign function 

sgn(ZV) = , DW = 75 sgn(H w ) , H w = 75 D w , (8) 

we write 

N 

sgn(H w ) = s g n (Vi)Xixl + P±sgn(H w ) . (9) 

i=l 

The first part is exactly handled in the subspace of eigenvectors x% °f %) the rest acting 



7l| satisfying 



in the orthogonal subspace is approximated by a minmax polynomial 

^minmax^) J=\ < 6 for p% +1 < X < /i^ ax , (10) 
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FIG. 2: The effective range function F(r) as a function of r/a on the 16 3 x 32 lattice at = 8.45 
for p = 1.4, together with an exponential fit. 

such that 

TV 

sgn(H w ) = s g n (^)xal + P± H w P minmax (if^) . (11) 

i 

The lowest iV = O(50) eigenvalues /ij of H w and the corresponding eigenvectors x% have 
to be computed. The degree of the polynomial P m i nmax depends on e and is proportional 
to the condition number k = /i^ ax / yU^ +1 of in the orthogonal subspace spanned by the 
modes {xi \ (1 — P±)Xi = 0}. For configurations generated with the Luscher-Weisz action, 
the condition number of Hw is essentially smaller in comparison with the Wilson action and 
only slightly depends on p. 

D(0) as described above is a massless Dirac operator. The massive overlap Dirac operator 



is easily obtained from the massless one, as 



D(m) = ( 1 - — ) D i°) + m ■ ( 12 ) 
Later on, in our spectral analysis we use the improved massless Neuberger operator 72, 
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£> imp (0) = (\ - D(0)j D(O) . (13) 

The resulting eigenvalues A imp of D imp (0) are the stereographical projection of the eigenvalues 
Ajy of D(0) from the Ginsparg- Wilson circle to the imaginary axis. The no n- zero eigenvalues 
appear in pairs A; mp = ±i A. In terms of these imaginary parts A, the spectral density and 
the characteristic properties of the eigenmodes are discussed. Also any mode cut-off A cut 
refers to this A. 



Finally, the improved massive overlap Dirac operator [43], which is used for example for 
the pion propagator in this paper, is 

D imp (m) = D imp (0) +m . (14) 



C. Eigenmodes of the Neuberger operator 

In principle, the emerging Neuberger overlap operator D(0) could have n_ plus n + exact 
zero modes, D(0)i/j n = 0, with n_ (n+) being the number of modes with negative (positive) 
chirality, jsipn = — 4>n (js^n — +VVi)- The index of D(0) agrees with the topological charge, 
Q = ri- — n + . In practice, however, there are always only negative or positive chirality 
eigenmodes for a given lattice configuration. The non-zero eigenvalues A^r and their non- 
chiral eigenvectors with D(0)ip\ N = Xn^Xm appear in complex conjugate pairs Ajv and X* N 
and satisfy J2 X (*P{ N ( x ),l^x N ( x )) = °- 

Locally, the eigenmodes ip\ of the Neuberger operator can be characterized by the scalar 
density 

Px (x) = \Mx)\ 2 = Y,^ C *( x ^ C ( x ^ Z>A(aO = l, (15) 

cr,c x 

where a and c denote spinor and color indices. 

Besides the scalar density, another important density needed is the pseudoscalar density 

Px 5 (a) = c * (^7f CT >A C W , ^Pas(x) = ±1 or (16) 

cr,a',c x 

for chiral or non-chiral modes, respectively. While the zero modes are globally chiral, the 
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non-chiral non-zero modes may still have a local chirality, 

Px± (x) = $>r "Vf c (x) , (17) 

CT,C 

with the projectors P± — (1 ± 75) /2 onto positive and negative chirality. 

D. Overlap definition of the all-scale and the mode-truncated topological density 



As for any 75-Hermitean Dirac operator satisfying the Ginsparg- Wilson re 
topological charge density for the Neuberger operator D(0) can be expressed as 



ation, the 



q{x) = -tr [ 75 (l-~D(0;x,x)^ , Q = ^g(x), (18) 



where the trace tr is taken over color and spinor indices. 



To compute the topological charge density, we use two different approaches [3l|, l33 |. 
In the first approach, the trace of the overlap operator is directly evaluated according to 
Eq. (j!8p . This is computationally very demanding and is therefore performed on only 53 (5) 
configurations in the case of a 12 3 x 24 lattice (16 3 x 32 lattice, respectively) at (5 = 8.45. 
The all-scale density q(x) computed in this way includes charge fluctuations at all scales. 
For this density, interesting anisotropic, global structures have been detected and discussed 



by Horvath et al. 34j . In this paper we will give more details on the multifractal properties 
of q(x). 

The second technique involves the computation of the topological charge density based 
only on the low-lying modes of the overlap Dirac operator. This approach is a gauge invari- 
ant filtering that leaves the lattice configuration unchanged (as well as the effective lattice 
spacing as the physical scale). Using the spectral representation of the Dirac operator, the 
truncated eigenmode expansion of the topological charge density reads 

Qx c Jx) = - (1 - ^) PasOO , (19) 

|A|<A cu t 

with p\s(x) defined in Eq. ffToT) . Truncating the expansion at A cut acts like an ultraviolet 
filter by removing the short- distance fluctuations from q(x). We will study in more detail 
how the properties of this density depend on the choice of A cut . With a suitable cut-off 
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the clusters of this density will coincide with the -R-clusters of (anti-)selfdual domains to be 
introduced in Section IV CI Note that the total topological charge Q = J2 X ^-W (• r ) * s no ^ 
affected by the choice of A cut , the level of truncation. 

E. Cluster analysis 

In this paper we will describe the properties of different densities in terms of their cluster 
properties. This is the appropriate place to describe the cluster algorithm in general terms. 
As an example we discuss this for the topological charge density q(x). 

1. The cluster analysis first requires tagging the lattice sites that will form the clusters. 
In the present case these are all sites with \q(x)\ > q cnt . 

2. Next these tagged sites are assigned to a set of link- connected clusters. Two sites x 
and y are said to be link-connected if they are immediate neighbors on the lattice 
connected by a link, x — y ± /t. Two link-connected and tagged sites x and y then 
belong to the same cluster unless q(x) and q(y) have opposite sign. Obviously, the 
latter veto is ineffective in the case of a positive density like the scalar density of 
eigenmodes. 

3. Cluster percolation is defined by the cluster correlation function which is given as the 
following ensemble average: 

fM _ ( E c e c (s) Q c (y) ) 6(r - \x - y\) 

where G c (x) is the characteristic function of a cluster c, i.e. Q c {x) — 1 if x £ c 
and c (x) = otherwise. In the definition (f20|) r is the Euclidean distance. If 
r max is the largest distance possible on the periodic lattice, we call the value /(r max ) 
"connectivity" . The onset of percolation is defined as the value of the cluster-defining 
quantity, i.e. of the cut-off g cut = g pe rc, such that the connectivity /(r max ) ^ for 

Qcut *^ *?perc- 

4. An important quantity describing the system of clusters is the Euclidean distance be- 
tween two clusters, for example the two biggest clusters. The distance d(c, c') between 
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two clusters labeled c and d can be defined as the maximum over sites x G c of the 
minimal Euclidean distance between x and any site y G d, 



5. Other quantities describing the system of clusters in a given configuration are the 
fraction of occupied volume, the fraction of the largest cluster to the total occupied 
volume, the size and charge distributions of clusters and the total multiplicity of 
clusters. 

We will use this terminology when we shall discuss qualitative features of certain densities. 
F. Random walkers 

For the characterization of a g-cluster with respect to its fractal dimension, we use the 
random walker method briefly described below. We consider a collection of independent 
random walkers starting from the maximum g max of the topological charge density inside 
a cluster among a collection of clusters defined by a cut-off q cnt . Depending on the ratio 
<?cut/<?max, the cluster is more or less extended. The cluster determines where the random 
walkers are permitted to go. The random walk is defined by the condition that at each 
time step r — > r + 1 the walker is required to jump with equal probability to one of the 
neighboring sites that also belong to the cluster. If g cu t /<?max is large, the cluster is small, 
such that the random walk soon arrives at a stationary regime. If (/cut/9max is smaller, the 
cluster is extended in one or more directions, the number of which one want to explore. 
The fractal dimension d* which is open for the random walkers will then be reflected by the 
return probability to the starting point, 



following a power-like decay with the number r of time-steps. This can be defined for one 
cluster or as average over all clusters. 

Also the distribution of all walkers in one cluster at a given time-step can be useful to 
characterize the shape of the cluster. This method is applicable to fractal clusters of any 
kind, for example also those of the scalar density of the individual fermionic modes. 
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FIG. 3: The Ginsparg- Wilson circle and the analyzed part of the unimproved and improved spec- 
trum for a Q = 3 configuration on the 16 3 x 32 lattice generated at f3 = 8.45. The insert shows 
the lowest part of the spectrum magnified. The zero eigenvalue is threefold degenerate. 

III. SPECTRAL RESULTS BASED ON THE LOWEST O(150) EIGENMODES 

A. The Ginsparg- Wilson circle and the topological charge from the index of D(0) 

Using the above construction of the Neuberger overlap operator, we perform its diago- 
nalization by a variant of the implicitly restarted Arnoldi algorithm. Mostly an amount of 
O(150) eigenvalues and eigenvectors per configuration has been computed and stored. The 
complex- valued eigenvalues A^r are located on the Ginsparg- Wilson circle of radius p around 
the point (p, 0) in the complex plane. This is shown in Fig. [3l 

In Fig. |4] (a) to (e) the distribution of topological charge Q = n_ — n + over the various 
ensembles is shown as determined from the number and chirality of zero modes of the lattice 
configurations. As has been said, we have never found zero modes with different chirality 
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simultaneously in the same configuration. 

From the distribution of topological charges shown in Figs. 0] the topological susceptibil- 
ities 



Xtop 




(23) 



have been evaluated. They are also given in Table I. It is remarkable that all physical lattice 
sizes realized with f3 = 8.45 are big enough to get the topological susceptibility equal within 
error bars. At the same time, the topological susceptibilities of the coarser lattices are only 
slightly larger, in accordance with the fact that other quantities are found to scale well with 
(5 [43j. The topological susceptibilities at (3 = 8.1 and (3 = 8.0 are equal to each other within 
error bars although the physical volumes differ by a factor five. All results are smaller than 
the topological susceptibility of quenched SU(3) gauge theory reported in Ref. [rjj to be 
(191 ± 5MeV) 4 . This value was given using fx = 160(2)MeV to set the scale. 



B. Spectral density of non-zero modes 

The spontaneous breaking of chiral symmetry by the dynamical creation of a nonvanishing 
chiral condensate is related to the spectral density p(A) of non-zero modes of the Dirac 



operator near zero by the Banks-Casher relation 75J] = — tt p(0). 

The spectral density of the continuous modes at finite volume is formally given by 

P(^V) = ^(Y1^-V)i ( 24 ) 

A 

where the sum extends only over positive (non-zero) values of A = Im Ai mp of the eigenvalues 
of the improved massless overlap operator D imp (0). The average (. . .) is taken over the 
ensemble of gauge field configurations. 

In the finite volume and for small eigenvalues the spectral density can be computed from 
the chiral low-energy effective theory. For A < Et, Et being the Thouless energy Ex ~ 
f^/YlyV, the low-energy effective partition function is dominated by the zero momentum 
modes, and the zero-momentum approximation of the chiral low-energy effective theory is 
equivalent to chiral random matrix theory. 

In random matrix theory the spectral density is given as 
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FIG. 4: The normalized distributions of topological charge Q, in subfigures ordered according to 
descending /?, for (a) 12 3 x 24 at (3 = 8.45, (b) 16 3 x 32 at [3 = 8.45, (c) 24 3 x 48 at (3 = 8.45, (d) 
12 3 x 24 at [3 = 8.10, (e) 16 3 x 32 at (3 = 8.00. Note that the physical volumes corresponding to 
(b) and (d) are roughly equal, and similarly for (c) and (e). 
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FIG. 5: The spectral densities a 3 p(X,V) at (3 = 8.45 for the three lattice sizes under study 
(12 3 x 24, 16 3 x 32 and 24 3 x 48) and the spectral densities at f3 = 8.1 and (3 = 8, for lattice sizes 
(12 3 x 24 and 16 3 x 32, resp.), together with fits using random matrix theory predictions. 



p(A, V) = S cfr W (Q) PQ&csVX) , (25) 
Q 

where S eff is an effective value of the chiral condensate, and 

Pq( x ) = 2^ J \Q\^ ~ ■ 7 IQI+i( a J IQ|-i( a 0) ' ( 26 ) 

is the microscopic spectral density [70] in the sector of fixed topological charge Q, expressed 
in terms of Bessel functions J n {x). 

We take the weights w(Q) of the sectors of topological charge Q (normalized to 
w(Q) = 1) from our measured charge distributions presented in Fig. HI To take into ac- 
count the effects of higher orders in chiral perturbation theory we also add a term aiA + c^A 2 
to the fitting formula ( 1251) . 

The spectral densities for the three different volumes available at the same j3 = 8.45, and 
the coarser ensembles of f3 = 8.10 (12 3 x 24) and f3 = 8.00 (16 3 x 32), together with our fits 
are presented in Fig. 5. 

In the region of small eigenvalues (A < 100 MeV) one can see a strong volume dependence 
of the spectral density. This is in agreement with predictions from quenched chiral pertur- 
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bation theory, where it has been shown 60J that the effective value of the chiral condensate 
S c g used in (125]) diverges logarithmically as the volume is sent to infinity. 



C. Localization and fractal dimension of the eigenmodes 

In this paper, we concentrate on details of localization and fractal dimension of the overlap 
eigenmodes. Later on, in Section ITVt similar aspects will be discussed for the topological 
density as derived from the overlap operator. 



A useful measure to quantify the localization of eigenmodes 44|, |47J] is the inverse parti- 
cipation ratio (IPR) 

/(A) = L 3 S L t J 2 (A) = Ll Lt 5>a(*) 2 , (27) 

X 

with a scalar density p\(x) that is normalized for all eigenfunctions to Y2 X P>^( X ) = 1- While 
one would have I = L\L t if the scalar density would have support only on one lattice point, 
the IPR decreases as the density becomes more delocalized, reaching 1 = 1 when the scalar 
density is maximally spread over all lattice sites. For independent Gaussian distributed 
ip a c (x) at each site (subject to an overall normalization) the density is not maximally 
spread, but still delocalized. This case corresponds to / = tt/2. 

The zero modes, that do not contribute to the chiral condensate, are always highly 
localized compared to the bulk of non-zero modes. But it turns out that the lowest non-zero 
modes may be localized as well. In Fig. [BJ in a series of subsequent timeslices, the scalar 
densities of the lowest (a) and of the highest analyzed non-zero mode (b) are shown for one 
Q = configuration 1 of the 16 3 x 32 ensemble generated at (3 = 8.45. The lowest non-zero 
mode is clearly localized in a few timeslices. But the scalar density of the highest analyzed 
mode is also slightly inhomogeneous. 

We study now more systematically the values of the IPR that occur in various parts of 
the spectrum. Also the zero modes are included in this discussion. In Fig. [7] histograms 
are shown with respect to the IPR for the complete set of all analyzed eigenmodes on the 
considered configurations for the five lattice ensembles, separately in bins of A. One can see 
again that the zero and lower modes are remarkably more localized, whereas the bulk of the 



1 This configuration serves as an example for many other structural observations in this paper. 
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(a) 



(b) 



FIG. 6: The scalar density of two typical non-zero eigenmodes of the overlap operator, shown in 
two-dimensional profile for a series of subsequent timeslices of a Q = configuration on the 16 3 x 32 
lattice generated at j3 = 8.45. (a) shows the first non-zero mode, (b) the highest (144th) analyzed 
mode. 

higher modes is delocalized. 

In Fig. [8] the IPR averaged over bins with a bin width AA = 50 MeV and for the zero 
modes (which are considered separately) is plotted. The average IPR shows a dependence 
on L s and on the lattice spacing a for zero modes and non-zero modes only in the range 
A < 150 MeV (for the first three bins). Beyond that interval the average IPR is practically 
independent of L s and a. There one finds (I) < 2.0 . It would be difficult to define an 
exact "mobility edge" through the localization. A proper definition would require to tell the 
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(a) 





(d) (e) 



FIG. 7: Normalized histograms of IPRs as functions of A, in subfigures ordered according to 
descending (3: (a) 12 3 x 24 at (3 = 8.45, (b) 16 3 x 32 at (3 = 8.45, (c) 24 3 x 48 at (3 = 8.45, (d) 
12 3 x 24 at (3 = 8.10, and (e) 16 3 x 32 at f3 = 8.00, for zero modes and non-zero modes in A bins 
with width 100 MeV. 
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FIG. 8: The average IPR for zero modes and for non-zero modes in A bins of width 50 MeV for 
the five ensembles. 

minimal number of modes above which the pion starts propagating. 2 

The volume and a dependence for the zero modes and the first three bins of non-zero 
modes allows some conclusions concerning the (fractal) dimension and the special localiza- 
tion properties of these modes. In view of a generalization of the IPR that will be made 
in a next step, we denote the effective dimension of an eigenmode based on the usual IPR 
as d*(2) (to remind us of the second moment of the scalar density p(x)). From the volume 



dependence of the average usual IPR, following (I) oc (V/a 4 ) 1_d * (2)/4 N, one gets this di- 
mension as an upper effective dimension. Similarly, the lattice spacing a dependence can be 
described by (I) = bo + b\ a d *( 2 ) -4 . From the first bin in Fig. [8] it is clearly seen that the IPR 
grows with the physical volume V at fixed a and with a -1 at fixed physical volume V. We 
fit the zero modes and the lowest non-zero modes in a couple of bins with a width AA = 50 
MeV. The results of the fit are presented in the uppermost curve of Fig. [9j For the zero 
modes an effective dimension of d*(2) = 2.2(1) and for the first interval of non-zero modes 
with A < 50 MeV an effective dimension d*(2) = 3.3(1) are found. 

Gubarev et al. 



46, 



47| came to somewhat different conclusions guided by the a dependence 



2 Although leading to a quantitatively wrong propagator, we will see that even the zero mode contribution 
alone would allow the pion to propagate. 
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of the usual IPR (/) = V I2 in pure SU(2) Yang-Mills theory, simulated with the Wilson 
action. Their fits for the a dependence of the IPR, with d*(2) restricted to integers, gave 
the lowest % 2 per degree of freedom for the choice d*(2) = 1 for the zero modes and choosing 
d*(2) = for the lowest non-zero modes. This discrepancy is most probably to blame to the 
Wilson action used to create the quenched configurations resulting in the proliferation of 
dislocations [52 ] . Our results concerning the usual IPR are similar to the findings of Aubin 



et al. 



44|. 



While the usual IPR of the low modes strongly depends on V and a, and hence the 
modes are localized with an effective dimension between d*(2) = 2 and 3, the IPR of the 
bulk of higher modes is independent of V and a. That means that these modes almost 
freely extend throughout d*(2) < d = 4 dimensions. We will see, studying the percolation 
behavior, that also the higher modes are not like simple plane waves, but experience a kind 
of soft localization. At this place it might be helpful to present in the Table II the average 
number of non-zero modes below some frequently used cut-off values A cut 3 . 



p 


lattice size 


a ^cut 


Acut 


average # of NZM 


8.45 


16 3 x 32 


0.1064 


200 MeV 


13.1(2) 






0.2128 


400 MeV 


36.2(2) 






0.3374 


634 MeV 


92.3(2) 


8.10 


12 3 x 24 


0.14 


200 MeV 


11.0(2) 






0.28 


400 MeV 


28.1(2) 






0.42 


600 MeV 


59.3(2) 






0.56 


800 MeV 


117.3(4) 



TABLE II: The average number of non-zero modes (NZM) below various spectral cut-offs for two 
ensembles used in this study. The cut-offs are given both in MeV and in 1/a. 

In the theory of the metal-insulator transition tools for a quantitative description of 

3 In the following the cut-off values 200, 400, 600 and 800 MeV on the 12 3 x 24 lattice at f3 = 8.1 should 
be considered as approximate values. Using the interpolated value for the lattice spacing at (5 = 8.1, the 
exact values are 195, 389, 584 and 778 MeV, respectively. 
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critical level statistics are given and the multifractal properties of the eigenfunctions are in 



73, 



78| to 
(28) 



the focus of interest. For this purpose, the notion of IPR has been generalized 

T m _ E,,J^(x)|^e £ (A-A n ) 

P[) En©e(A-A n ) 

with a free parameter p and a window function e (A — A n ) = 1/e for |A — A n | < e/2 and 
vanishing elsewhere. Depending on whether the metal phase, the insulator phase or the 
critical region is met, this quantity would scale in a different way with the volume of the 
specimen (see Table III, results see Fig. [9]). 



I P (X) oc L~ d( - p -V 


metallic 


electrons propagate freely 


I P (X) oc L^pXp- 1 ) 


critical 


electrons propagate along 
low-dimensional structures 


ip(A) = const 


insulator 


electrons do not propagate 



TABLE III: Volume scaling of generalized IPR's according to the metal-insulator transition analogy. 



For the "metallic" phase, d is the embedding dimension of the system, in our case d = 
4. For the critical region the effective dimension becomes d*(p) < d. In this sense, all 
eigenstates with A < 200 MeV definitely belong to a "critical region". From the standard 
IPR I = L t I 2 (p = 2) we have concluded that the dimension of the zero modes is close 
to two. In the next bins of 50 MeV each, the dimension quickly exceeds three, before it is 
finally approaching four. 

For the "critical region" and up to A = 400 MeV we present now the results for the higher 
order p fractal dimensions d*(p), which allows to recognize a varying type of localization at 
various levels (heights) of the scalar density p(x). Notice the decrease of the multifractal 
dimension with increasing p, resulting from the fact that I p explores regions of higher scalar 
density. For p = 4, for example, one sees generally a reduced dimension compared to p = 2, 
interpolating between "filamentary" (d*(4) « 1 for zero modes) and "surface-like" (d*(4) « 2 
for the next two bins with A < 100 MeV), becoming <i*(4) rs 3.5 at higher A. This indicates 
that the regions of higher scalar density are geometrically distinct from those regions where 
an eigenmode is present only with a tail of the scalar density. For p = 20 - 30 the analyzed 
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FIG. 9: The fractal dimension d*(p) obtained from fits of the volume dependence of the averages 
of the generalized I p , presented for zero modes and for non-zero modes in A bins of width 50 MeV 
for the three ensembles with different volume and common (3 = 8.45 . 

fractal dimension does not grow anymore. The envelope tells that, with respect to the 
maxima of the scalar density, zero modes and non-zero modes up to A = 100 MeV are all 
characterized by a less than one- dimensional localization [i.e. isolated peaks), whereas in 
the following spectral region two-dimensional or three-dimensional localization prevails for 
the regions of highest scalar density (cf. Fig. [6]). In contrast to this, the standard IPR based 
on J 2 alone cannot resolve these details and would describe the modes in the spectral region 
A > 100 MeV as four-dimensional. 

It is tempting to conjecture a pinning-down of the low-dimensional low-lying modes to 
specific confining objects (vortices, monopoles, close-by meron pairs etc.) and a relation 
to the localization of topological charge (see Section HVl in this context). It is important 
to recognize that the difference between the zero modes and the lowest (say, 10) non-zero 
modes is less pronounced than it seems in the result of binning. In marked contrast to this, 
the highest analyzed non-zero modes are really qualitatively different. This could already 
be concluded from Fig. El 

We show in Fig. [TCD (a) the result of a cluster analysis 4 of 170 configurations on the 16 3 x 32 
4 The cluster analysis with respect to some observable was introduced in Section III El 
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FIG. 10: Cluster analysis of some individual eigenmodes as listed in the legend box, averaged over 
170 configurations on 16 3 x 32 generated at /3 = 8.45. The Pcut/Pmax dependence is shown of (a) 
the total number of separate clusters at the height p C ut, (b) the connectivity as a signal for the 
percolation of the mode, (c) the missing norm of the mode and (d) the effective dimension of the 
mode at height p cu t determined by the return probability (|22p of random walkers. 

lattice generated at (3 = 8.45. The average is over the zero modes in this subensemble, the 
1st, the 10th, the 30th etc. up to the (arbitrarily chosen) 120th non-zero mode. We analyze 
the average cluster composition as function of the cut-off scalar density p C ut/Pmax (relative to 
the maximum of the scalar density of the respective mode p max = ™-&x x p\(x)). In this type 
of cluster analysis one attempts to decompose the set of lattice points with p\(x) > p cut into 
one or more connected clusters, each enclosing secondary maxima of p\(x). For the highest 
analyzed mode, beginning at p cut = 0.9 p max , rapidly further maxima become visible, and 
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the number of separate clusters grows to ~ 20 at p cut = 0.35 p mSLX . In contrast to this, the 
typical zero mode only slowly develops up to ~ 5 maxima at p cnt = 0.1 p max . From the zero 
mode to the highest analyzed mode there seem to exist several centers that "attract" the 
mode, but the mechanism might be different. The lowest modes might be pinned down to 
some lumps of topological charge (o la Diakonov and Petrov 7^|, whereas the localization 
of the highest (analyzed) modes might be the result of Anderson-like localization in the 
random gauge field background. 

The lowest ten non-zero modes are very similar to the typical zero mode. The zero modes 
percolate, i.e. extend over the full lattice, only at a height below p cut = 0.25 p max - 

The highest analyzed non-zero mode percolates at p cut = 0.4 p max (with the maximal 
scalar density being smaller, of course). This can be seen in Fig. [10] (b). The difference in 



the percolation behavior can be concluded also from the missing norm 80] 



{a;|p(x)>p cut } 

as function of p cut /p max that is shown in Fig. [10] (c). Here the similarity between the zero 
modes and the first 10 non-zero modes is similarly clear. 

At each height level of p\ the effective shape of the leading cluster of the respective mode 
can be explored by the random walk method described in section Hi Fl In the present context 
one studies how the return probability of the walkers to the maximum of p\(x) decreases 
with an increasing number r of steps. A fit of the power law decrease ([22]) provides the 
effective dimensions d* shown in Fig. [TU](d). We notice that the average dimension averaged 
over the 10 lowest modes continuously rises from d* = to d* = 3.5. It becomes three- 
dimensional at the same p C ut/Pmax when the missing norm begins the final steep drop from 
80 % to %. As for the 120th non-zero mode, the effective dimension d* rises more steeply 
and reaches d* ~ 3 at higher p C ut/Pmax when the missing norm is already dropped to 60 
% . For a survey of dimensionalities of overlap eigenmodes see Tables VII and VIII in the 
Appendix. 
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D. Chiral properties of the non-zero modes 



In this section we characterize the different parts of the low-lying spectrum of non-zero 
modes by the distribution of sites over the local chirality. Later on, the local chirality of the 
lowest modes is shown to be locally correlated to the topological density 

In a given lattice point x, one might be interested in the intensity ratio 

r \ 0) = 30 

V\-[x) 

of the two chiral projections for a given mode ip\. For the chiral zero modes this ratio is oo 
or everywhere. For the non-chiral, non-zero modes it is some function of x. It is useful to 

n 

relate r\(x) to [50| 

X x (x) = - arctan (v^aM) - 1 e [-1, +1] . (31) 

It will be seen that relatively low-lying modes have regions where the local chirality carries 
information about the (anti-)selfdual character of the background field, to be precise, where 
0.5 < \X(x)\ < 1. 

In Fig. [TT] and Fig. [12] normalized histograms with respect to the local chirality X(x) 
are shown for the Q = subsample (containing 37 configurations) of the 16 3 x 32 lattices 
generated at (3 — 8.45 with the lattice sites selected according to the scalar density p(x). The 
left figures show the 10 different X-histograms for the lowest 20 individual non-zero modes 
(the pair with A and —A gives rise to the same histogram), whereas the right figures show 
the histograms where all non-zero modes belonging to a A bin (with a bin size 100 MeV) 
are contributing to 8 histograms covering the spectral range up to 800 MeV. The different 
rows show different cuts applied restricting the lattice sites x to those with p\(x) > p cut . 
The cuts are chosen such to represent 1 %, 6.25 %, 12.5 % and 50 % of all lattice sites. 

For zero modes (not included here) the histograms would be (^-functions at X = —1 or 
X = 1 independent of any cut with respect to the scalar density. Although being non-chiral 
when integrated over space-time, the first non-zero modes locally may have a considerable 
degree of chirality (peaking in X somewhere between 0.75 and 0.9) for the leading 1 % or 
6 % of the lattice sites with the highest scalar density. It turns out that with a higher, more 
restrictive cut the degree of local chirality can be enhanced. A similarly stronger effect can 
be seen towards the lower modes. That means that the 5-10 maxima of the scalar density 
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FIG. 11: Normalized histograms of the local chirality in the Q = subsample (consisting of 37 
configurations) of the 16 3 x 32 lattices generated at j3 = 8.45. Left: for the lowest ten pairs (with 
positive and negative A) averaged over the subsample. Right: for all modes from the subsample, 
averaged over A bins of width 100 MeV. Upper row: with a cut for 1 % of the sites with biggest 
scalar density p(x), lower row: with a cut for 6.25 % of the sites with biggest p(x). 
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(c) (d) 

FIG. 12: The same as Fig. [TT] but with less exclusive cuts with respect to the scalar density: Upper 
row: with a cut for 12.5 /lower row: with a cut for 50 % of the sites with biggest p(x). 
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of the lowest non-zero modes are localized in space-time region where the local chirality is 
strongly correlated to the topological charge density of the gluonic background. The highest 
modes, as we have found, have O(20) local maxima. In contrast to the lower modes, the 
histogram of local chirality for the higher ones is dominated by a broad bump of low chirality 
X(x) distributed around zero for all cuts in p(x). This means that the maxima (localization) 
of the highest modes are not the result of pinning down at spots of topological charge density, 
rather they are the result of waves scattering on a random gauge field background. 

One can see that the eigenmodes with eigenvalues up to 200 MeV show an enhancement 
of local chirality around X = ± (0.5. ..0.75) as long as one focusses on the 6 % of the lattice 
points with the strongest scalar density. If 50 % of the lattice sites are included in the 
analysis, only the first two pairs of non-zero eigenmodes still show some enhancement of 
local chirality around X = ±0.75. 

We will later (at the end of Section IIVEI) consider the correlation function between the 



topological charge in the gluonic background and the local chirality X(x) of individual 
non-zero modes. The result establishes that they are interrelated over some distance in 
space-time. 

IV. TOPOLOGICAL CHARGE DENSITY FROM THE OVERLAP OPERATOR 

A. Two-point function of the all-scale and the mode-truncated density 

The topological charge density of the gluonic field seems to play a deciding role for the 
localization of the lowest modes. In this Section we are studying the properties of the 
overlap definition of the topological density according to Eqs. (Tl8|) and (TT9"]) . The striking 
difference between the all-scale topological charge density and the mode-truncated one can 
be illuminated studying the point-to-point correlator of the topological density. Formally, 
the definition of the two-point function is the same for q(x) and q\ cut (x): 

C qq {r) = i=-_ — - < for r > . 32 

22x,y S ( r - \ X -V\) 

Reflection positivity, i. e. the pseudoscalar nature of q(x) together with positivity of the 
metric in Hilbert space, demands that the topological charge density correlator is negative for 
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FIG. 13: The topological charge density correlator for the 12 3 x24 lattice configurations at (3 = 8.10. 
(a) The correlator of the mode-truncated density q\ cut (x) for various A cut compared with the 
correlator of the all-scale density q(x), measured only on 53 configurations. The right figure (b) 
magnifies the region in r where the mode-truncated correlators become negative for sufficiently 
large A cut . 

all distances, C qq {r) < for r > 0. Of course, the vacuum expectation value ( q(x) q(y) ) 
must be a function only of the distance r = \x — y\ > 0, due to Euclidean rotational 
invariance. 

While the topological charge correlator C qq (r) in the continuum should be negative for 
nonvanishing distance, the topological susceptibility, which is the integral over the correlator 
Xt — J dx ( q(0) q(x) ), must be positive. To solve the apparent contradiction, formally 
divergent contact counterterms of the form C qq (x) — » C qq (x) + c\8{x) + C2A5(x) + csA 2 5(x) 



have to be introduced in the continuum theory 

Since overlap fermions are not ultralocal, the fermion action cannot be strictly reflection 
positive. Because the overlap operator D(0) has a finite range in space-time (see Fig. [2]) 
practically independent of /3, we expect that the correlator C qq (r) will have a positive core 
with a fixed width in lattice units. The negative tail beyond the core radius corresponds to 
the above mentioned theoretical requirement. This has been discussed recently in Ref. 

For the correlator of the mode-truncated density q\ cut (x) one would expect that this 
structure is smoothed out and that below some A cut the corresponding two-point function 
will be positive at all distances. In Fig. [13] we show the correlator of the mode-truncated 
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FIG. 14: Comparison of the topological charge density correlator based on the all-scale density 
q(x) for the 12 3 x 24 lattice at (3 = 8.10 (53 configurations) with the same correlator for the 16 3 x 32 
lattice at (3 = 8.45 (only 5 configurations) and (5 = 8.60 (2 configurations). The right figure (b) 
magnifies the region in r where the correlator becomes stronger negative when one gets closer to 
the continuum limit. 

topological charge density for various A cut for the 12 3 x 24 lattice configurations generated 
at (3 = 8.10. In Fig. [13] (a) they are compared with the correlation function of the all-scale 
topological density measured on 53 configurations. In Fig. [13] (b) we show the correlator 
only for different mode-truncated densities, however magnified, in the region of r where they 
start to develop a negative tail for sufficiently high cut-off. This happens at R > 6 a for 
A cut = 800 MeV, at R > 7 a for A cut = 600 MeV, and at R > 9 a for A cut = 400 MeV. 
As long as it were positive, this correlator would be interpreted in terms of the number 
density and the convoluted profile of effective topological charge clusters. It is interesting, 



that higher resolution in the present case (similarly to doing less smearing 82j or doing 



restricted cooling [2J, [83j) results in effectively attractive correlations between opposite 
sign topological charge clusters not only in full QCD but also in quenched gauge fields as 
emphasized in Ref. |82j. Since the total topological charge does not depend on A cut , the 
growing negative tail must be compensated in the susceptibility ^2 x C qq (\x\) by a growing 
positive core with increasing A cut . 

In order to demonstrate how the correlation function of the all-scale topological charge 
density varies with the lattice spacing we compare in Fig. [TH the correlator for 12 3 x 24 
lattices at /3 = 8.10 (53 configurations) with the correlator measured on 16 3 x 32 lattices 



at (3 — 8.45 (only 5 configurations) [6|] and 16 x 32 lattices at j3 — 8.6 (2 configurations) 
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An extrapolation of the function a((3) to the highest (3 gives a(/3 = 8.6) = 0.096(1) . The 
volumes of the first two lattices are approximately equal. The required negativity of the 
correlator develops only for a sufficiently fine lattice. The much bigger number of clusters 
at higher values \q(x)\ ~ 0.25 g max (see Section HVCp seems to be necessary to achieve this. 
Thus, the cluster multiplicity of the all-scale topological density simply reflects the size of 
the lattice in lattice units. 

B. Two-dimensional profile seen with varying resolution 

As a vacuum expectation value, the point-to-point correlation function obviously has to 
be rotational invariant. Individual lattice configurations, however, are richer in structure 
and necessarily have a locally anisotropic structure. This results from lower- dimensional 
structures present in the configurations that locally break rotational invariance. 

As an example, a snapshot with four levels of resolution is presented in Fig. [15]of the same 
two-dimensional section through the typical configuration, already considered in Figs. [6] and 
[2H of the ensemble of 16 3 x 32 lattices at (3 — 8.45. The respective resolution can be inferred 
from the corresponding A cut . The two-dimensional section is periodically doubled in the two 
directions such that the one-dimensional structure (apparently one "ridge" per elementary 
cell in this two-dimensional world) can be recognized. In particular, this ridge is clearly 
visible at A cut = 200 and 400 MeV, being already less clearly visible at A cut = 634 MeV. 
At a resolution of 400 MeV and higher, also single peaks (a 0-dimensional structure in this 
two-dimensional world) on top of the ridges become visible. Describing the landscape from 
the perspective of the peaks of density in Fig. [T5] (c), the next-higher dimensional struc- 
tures ("ridges") are features becoming discernible only at lower height of the density. In 
a cluster analysis, with g cut chosen too large, one might miss the ridge catching only the 
peaks. Varying q cut the dimension revealed by the cluster analysis will be very reminiscent 
of the multifractal dimensions of the individual fermionic modes. Compared to these visible 
features of the three different mode-truncated densities, the cut through the all-scale topo- 
logical density is seemingly completely random. Nonetheless, structure can be discovered in 
the "full" density, too, simultaneously taking all four dimensions into account through the 
cluster algorithm, as we shall see next. 
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(a) (b) 




(c) (d) 

FIG. 15: Two-dimensional profile through the typical Q = configuration generated at (5 = 8.45 
on a 16 3 x 32 lattice. The profile is periodically doubled in both directions for better visibility of 
extended structures. The cuts show, from (a) to (c), the mode-truncated density for descending 
levels of truncation A cut = 200 MeV, A cut = 400 MeV and A cu t = 634 MeV. The all-scale topological 
density is shown in (d). Note the 20- fold larger vertical scale in (d). 

C. Cluster analysis of the all-scale topological charge density 

In order to make the structural analysis more quantitative, the cluster analysis as de- 
scribed in Section III El is performed. Let us begin with the all-scale topological charge 
density. The structure is changing with the threshold value q cut . As the reference value 
for g cut serves the maximum of the topological density in each configuration, g max . Thus, 
we take < g C ut/<?max < 1 as the running parameter. We remark that g max for the mode- 
truncated density is weakly changing with A cut (around 0.005 in lattice units). For the 
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A C ut 


200 MeV 


A r\r\ H If T T 

400 MeV 


600 MeV 


634 MeV 


800 MeV 


all-scales 


(3 = 8.10 
12 3 x 24 


5.04(14) 10~ 3 


6.45(14) 10" 3 


7.66(15) 10~ 3 




8.88(14) 10~ 3 


2.82(07) 10~ 2 


P = 8.45 
16 3 x 32 


3.51(17) 10~ 3 


4.45(18) 10~ 3 




5.39(19) 10~ 3 




1.97(06) 10~ 2 



TABLE IV: The ensemble averages of q ma , x a 4 for the mode-truncated topological charge density 
(for various A cu t) and for the all-scale density for two lattice ensembles. 

all-scale topological density g max is a few times larger (see Table IV). 

In Fig. [16] the gcut/<?max dependence of the cluster multiplicity, of the relative size of the 
largest cluster with respect to all clusters, of the distance between the two biggest clusters 
and of the connectivity are presented. This is shown for two ensembles of practically the 
same volume, 12 3 x 24 at (3 = 8.10 and 16 3 x 32 at [3 = 8.45, and, in order to study the 
(3 dependence alone, for a few configurations of the same size 16 3 x 32 in lattice units, but 
with a finer lattice spacing at (3 = 8.60. 

The most striking difference between the first two cases with almost equal physical volume 
is the roughly four times larger multiplicity of clusters that is reached on the finer lattice. 
In both lattices percolation (nonvanishing connectivity /(r max )) sets in close to where the 
multiplicity reaches its maximum. The percolation threshold g perc is practically independent 
of the lattice spacing. The percolation finally ends with g cut — > with only two, globally 
extended clusters of opposite sign of topological charge. With decreasing gcut/<?max, but still 
above the percolation threshold at g cu t/<?max = 0.2, the number of clusters increases strongly 
with decreasing cut, in particular on the fine lattice. Since number and size of all clusters 
increase, the relative size of the biggest cluster decreases. This trend is reversed beyond 
the percolation threshold where the relative size of the largest cluster rapidly reaches 50 %, 
practically independent of the lattice spacing. The distance between the two biggest clusters 
shown in Fig. [16] (c) is given in lattice spacings. Before percolation the distance is therefore 
only seemingly bigger on the finer lattice. Taking the ratio a{(3 = 8.10)/a(/3 = 8.45) m 1.4 of 
the lattice spacings (see Table I) into account, the average distance between the two leading 
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FIG. 16: Cluster analysis of the all-scale topological density: The q cut dependence is shown of (a) 
the total number of separate clusters, (b) the size of the largest cluster relative to all clusters, (c) 
the distance between the two largest clusters in lattice units and (d) the connectivity (see section 
III Cp . The data is plotted both for the 12 3 x 24 lattice configurations at f3 = 8.10 (53 configurations, 
in red), for the 16 3 x 32 lattice configurations at (3 = 8.45 (only 5 configurations, in green) and for 
the 16 3 x 32 lattice configurations at f3 = 8.60 (only 2 configurations, in blue) being averaged over. 

clusters is approximately the same and independent of qcut/qm&x all over the non-percolating 
regime. The finer the lattice, the faster the minimal distance between the two leading clusters 
is approaching the minimal distance of two lattice spacings as soon as g cut < q pe rc- In the 
final state the two clusters are close to each other everywhere, kept as separate clusters only 
because of the different sign of q(x). In other words, the remaining two global clusters of 
opposite charge are everywhere thin of 0(a) and entangling each other. The finite distance 
in lattice spacings is independent of (3 reflecting nothing else than the non-ultralocality of 
the overlap operator. 
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Cluster Nr. i Cluster Nr. i 

(c) (d) 

FIG. 17: Cluster analysis of the all-scale topological density: The fractional volume ((a)+(b)) 
and the total charge ((c)+(d)) of the corresponding clusters, ranked according to their occupied 
volume, is shown for two values of qcut / Qmax' in the percolated regime (left) for (?cut/<7max = 0.10 
and in the cluster-separated regime (right) for q cut / 'g max = 0.30. The plot represents one single 
configuration of the 16 3 x 32 lattice ensemble generated at j3 = 8.45. 

For the two lattices 16 3 x 32 with different lattice spacing we find a surprising similarity 
in all four plots. Even the distance between the biggest clusters coincides within large errors 
due to the low statistics (5 and 2 configurations for (3 = 8.45 and (3 = 8.60, respectively). 
In particular, the cluster multiplicity seems to depend mostly on the lattice size in lattice 
units. 

How the whole g-cluster composition with respect to charge and volume varies between 
the percolating regime (g cut = 0.1 g max ) and the non-percolating regime (g cut = 0.3 g ma x), 
is illustrated in Fig. [T71 As long as percolation has not set in, there is a broad spectrum 
of cluster volumes. The spectrum of cluster topological charges is surprisingly similar, up 
to the random sign. The ranking of clusters is according to the cluster volume, and the 
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FIG. 18: The probability of return to the cluster center x c for random walks, restricted to topologi- 
cal clusters of the all-scale topological charge density, shown as a function of the step number. The 
clusters are distinguished by different cuts with respect to the density that they enclose. The cases 
Qcut/Qmax = 0.4, 0.3, 0.2, 0.1, and 0.0 are considered on the coarser lattice 12 2 x 24 at (3 = 8.10. 
The results of the power fits (see the curves) are labeled by the effective dimensions d* . 

modulus of the cluster charges almost follows this ranking. Thus, the average charge density 
in each cluster is approximately the same. When percolation is completed, there are only 



two, oppositely charged clusters left over with approximately equa 
situation has been first discovered and discussed by Horvath et al. [33] 



volume. This peculiar 



D. Dimensionality and multifractality of the all-scale topological charge 

In this subsection a more explicit description of the fractal dimension and the multifractal 
properties of this topological density will be given. These properties are difficult to visualize 
in three dimensions because of the close packing and the huge multiplicity (at the percolation 
threshold) of the corresponding g-clusters. 

For our purpose we use the random walker method described in section III Ft In Fig. [18] 
we show in a double-logarithmic plot the decay of the return probability on the 12 3 x 24 
lattice. In an open d*-dimensional space the return probability after r steps would follow 
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the power law P(0, r) = l/(27rr) d */ 2 . It turns out, that by extending the clusters (i.e. by 
chosing (/cut /(/max lower and lower) the power of decay becomes larger. The corresponding 
effective fractal dimensions are ranging from d* = 2.7 for low (/cut/oWx ~ 0, over d* = 1.5 
at <7cut/<?max = 0.1 and d* = 1.1 at g cu t/<?max = 0.2 where the g-cluster multiplicity reaches 
its maximum, to d* = 0.7 at g C ut/<?max = 0.3 where the clusters are really restricted to the 
immediate neighborhood of the maxima. For a survey of dimensionalities of clusters of the 
topological charge see Tables V and VI in the Appendix. Even if the values of the fractal 
dimensions d* provided by the random walker algorithm do not exactly coincide with the 
dimensions obtained from other methods, the tendency of the results of the random walk 
method indicates that the landscape of the all-scale topological charge is multifractal and 
becomes nearly three-dimensional only far below the percolation threshold. Percolation sets 
in at g cut ~ 0.25g max , when the effective dimension is still low as d* = 1 — 1.5 . 

A second argument comes from a consideration of the cumulative cluster charge. One 
considers a collection of g-clusters primarily defined for a given cut-off g C ut/<?max- One covers 
the set of clusters with a four- dimensional sphere of radius R centered at the center x c of 
one particular cluster (i.e. the maximum of \q(x)\ inside the cluster c). One sums that 
part of the cluster charge that is located in the intersection of the cluster and the covering 
sphere. The relative cumulated cluster charge is the ratio between the cumulated charge 
to the full charge of the cluster and it depends on the radius R. The dependence of the 
relative cumulated cluster charge on the radius R will start with some power and finally 
reach unity at saturation. The observed power behavior can give information about the 
fractal dimension d* < d of a g-cluster. The initial growth with R may be parametrized 
as R d * , where d* describes the local dimension near the cluster center and depends on the 
cut-off gcut/<?max- One can see in Fig. [19] that a linear fit works very well as long as qwt/gmax 
is sufficiently high such that the clusters are isolated and do not percolate. The beginning 
percolation is characterized by an initial growth fitted like Q C umui ~ R 3 - 

One might expect problems when R is so large that the sphere covers many different three- 
dimensional branches of a g-cluster which might appear disconnected inside the sphere but 
are actually globally connected (and therefore counted in the cumulative charge according 
to the cluster membership). This would happen at very low cut-off q cnt when only the two 
largest clusters have survived. Only then the cumulated charge of each of the clusters would 
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(b) 



FIG. 19: The covering sphere method (see text): (a) for the coarse lattice 12 3 x 24 at = 8.10, (b) 
for the fine lattice 16 3 x 32 at (5 = 8.45. The fitted power behavior of the accumulated charge jumps 
from one-dimensional to three-dimensional at some percolation threshold (/ C ut/<?max (see text). 

grow like ^cumulative ~ R A although the local dimension is d* = 3 5 . 

One can see in Fig. [TH] that the transition from linear to cubic behavior happens for the 
coarser lattice (left panel) somewhere between g C ut/<?max = 0.1 and 0.2 . The actual power 
behavior for the lowest shown cut-off is like Q cumulative ~ R 3 ' 1 . On the finer lattice, for the 
higher cut-off g C ut/<?max = 0.2 the situation is already close to the d* — 3 regime. The fitted 
power behavior is like (^cumulative ~ R 25 - For g C ut/<?max = 0.1 it is like R 3 3 . 

For a survey of dimensionalities of the topological density (as well as of overlap eigen- 
modes) inferred from the covering sphere method see Table VI (and Table VIII) in the 
Appendix. 



E. Cluster analysis of the mode-truncated topological charge density 

Let us now discuss the effect of mode-truncation on the cluster structure of the topological 
charge density. We compare in Fig. [20] for the two ensembles, 12 3 x 24 at (3 = 8.10 (left) and 
16 3 x 32 at (3 = 8.45 (right) with practically equal physical volume, the ensemble averages 
of the cluster multiplicity, of the fractional volume filled by all clusters (packing fraction) 

5 This consideration was the motivation to propose the random walker method described in Section III Fl 
The walkers are not allowed to penetrate to or return from another branch of the same cluster. 
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(e) (f) 

FIG. 20: Comparison of the cluster structure of the mode-truncated density with three resp. four 
cut-offs A cu t with the all-scale topological charge density. The effect of (5 at the same physical 
volume can be seen by comparing the left side (12 3 x 24 lattices at j3 = 8.10) with the right side 
(16 3 x 32 lattices at j3 = 8.45). Top row: the q cu t dependence of the total number of clusters ((a) 
and (b)). Note the different multiplicity scale on the right for clusters of the all-scale topological 
charge density. Middle row: the packing fraction of all clusters together ((c) and (d)). Bottom 
row: the connectivity ((e) and (f)). 
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FIG. 21: Comparison of the cluster structure of the all-scale topological charge density with the 
mode-truncated density for three values of the cut-offs A cut : (a) the fractional size of the largest 
cluster and (b) the distance between the two largest clusters in lattice units, both as function of 
9cut/?max- The data refers to the 16 3 x 32 lattice at /? = 8.45 only. 

and of the connectivity /(r max ). In each panel we show this for 3 or 4 cut-offs A cut together 
with the same quantity denned for the all-scale density. Notice the different scale for the 
cluster multiplicity for the all-scale density on the right of the upper panels. For the coarser 
lattice presented on the left and the finer lattice presented on the right, the tendency is the 
same. The maximal number of clusters is much smaller for the mode-truncated density than 
for the all-scale density. It rises with the cut-off A cut . For a fixed cut-off A cut in MeV, it 
is approximately 40 to 100 % larger on the finer lattice than on the coarser lattice. Thus, 
the cluster multiplicity defined with a fixed ultraviolet cut-off A cut , is strongly discretization 
dependent. The maximum of cluster multiplicity is reached at ever smaller q cn t/qma.x for the 
finer lattice. The packing fraction, i.e. the fraction of volume filled by the clusters, for the 
all-scale density rises faster with decreasing g cu t on the finer lattice than on the coarser. On 
the other side, for the mode-truncated densities, the final rise of the packing fraction happens 
almost simultaneously for g cut < 0.1 g max for all cut-offs A cut under consideration. This means 
that, as long as g cu t/<7max > 0.1, the clusters remain well separated. These differences are 
washed out on the coarser lattice. On the finer lattice the maximal cluster multiplicity (or 
the percolation threshold) for the all-scale density is found at much higher g cu t/g m ax value 
than for the mode-truncated density including the highest investigated A cu t = 634 MeV. On 
the coarser lattice this happens to the all-scale and the mode-truncated density at the same 

max • 
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For the finer lattice only, in Fig. [21] the relative size of the largest cluster and the distance 
between the two largest clusters is plotted in dependence on g cu t/<?max- Each panel shows 
this for various A cut and for the all-scale density. The behavior for the different cut-offs A cut 
of mode-truncated densities is similar. The main difference of the mode-truncated density 
compared to the non-truncated (all-scale) density is that the relative size of the largest 
cluster remains relatively constant before and across the percolation threshold although the 
number and size of all clusters increases with decreasing g cut before percolation sets in. The 
largest cluster is always much bigger than the other clusters, and this effect is the stronger 
the lower the cut-off A cut is. It can be clearly seen that the distance d(a, b) between the two 
largest clusters a and b of the mode-truncated density does not drop rapidly to 2 lattice 
spacings immediately after the onset of percolation. This is because these clusters remain 
well-separated even at g cu t/(7max < 0.1. 

All these observations point towards the conclusion that the mode-truncated density 
allows to define (/-clusters of topological charge that, although they also begin to percolate 
at some (relatively low) height (cut-off q cn t), remain well-separated even at this level of 
the density. One of the clusters is much larger than the remaining ones. These properties 
they have in common with the i?-clusters that will be defined on the basis of the local 
(anti-)selfduality R(x) in Section IV CI Thus, these clusters allow an interpretation close 
to the traditional "lumpy" instanton or caloron gas picture, however with the tendency to 
coalesce into larger objects of same-sign topological charge density ("instanton clumping"). 
At higher A cu t an attractive correlation between opposite charge clusters appears. 

Given this cluster interpretation of the mode-truncated topological density, we return to 
the question whether the chirality X n (x) of the n-th non-zero mode is locally correlated with 
the topological charge density. To answer this question we consider the correlation function 
between the density gA cut iy) and the local chirality X{y) over some distance r, C q x{r; A cut , n) 

defined as _^ 

n ( , v ExJ 9A cut (^) x n(y) ) S(r- \x-y\) 

C qX {r; A cut ,n) = = — , ^ . (33) 

The result is shown in Fig. [22] for A cut = 100 MeV. The local chirality that enters the 
correlator in the figure is evaluated for the 2nd, the 6th, the 10th, the 30th, the 50th etc. 
up to the 120th non-zero mode. We keep in mind that this mode-restricted topological 
charge density is influenced (on average) by the lowest 9 or 10 pairs of non-zero modes. 
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FIG. 22: The correlation function between the mode-truncated topological charge density for 
^cut = 100 MeV and the local chirality of selected non-zero modes from the lowest one to the 120th 
mode. The correlation function is an average over the 37 configurations with Q = on the 16 3 x 32 
lattice generated at j3 = 8.45. 

We see that the positive correlation of the (collective) topological charge density on one 
hand with the local chirality of individual overlap modes up to the 10th mode on the other 
hand is strong and ranges over up to ten lattice spacings (similar to the distance where 
the correlation function C qq of gA cut changes the sign. At larger distance the local chirality 
gets anticorrelated to the mode-truncated topological charge density. Although weaker and 
over shorter distance, the 30th and higher modes are still positively correlated with the sign 
of the truncated topological charge density. For these higher modes the anticorrelation is 
absent at larger distance. 



V. THE INFRARED FIELD STRENGTH TENSOR 



A. Basic idea 

We have considered in Section [TV] the topological charge density as represented by the 
overlap operator according to Eqs. ( jl"8l) or ( fl"9l) . We have seen that the structure of the mode- 
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truncated topological charge density strongly depends on the number of included lowest 
eigenmodes. Restricting to fewer modes is tantamount to an ultraviolet filtering applied to 
the topological density. Our point of view as explained in Section [I] is that the low-lying 
modes represent the physically relevant degrees of freedom for hadronic physics. Thus, 
filtering should be suitable to highlight the important structures. Here we briefly discuss 
another filtering method similar to what was proposed by Gattringer and is suggested 
by the newly discussed representation of the field strength tensor through the overlap Dirac 
operator (39), |4(]] . The main purpose is to get an ultraviolet filtered, infrared field strength 
tensor. This would allow to assess the local degree of (anti-)selfduality. The local selfduality 
or antiselfduality of a gauge field is an important feature in all semiclassically motivated 
vacuum models, and the search for it in generic lattice configurations may play a similar (or 
superior) role compared to the search for lumps of topological charge, which was the main 
strategy of access before. 

The starting point is the observation that in the continuum the square of the Dirac 
operator D, projected with the help of a combination of 7^ matrices, represents the field 
strength tensor 

F^(x) = ~ tr Dirac (<v D 2 {x)) . (34) 



Here trpj rar is the trace taken over spinor indices. In Ref. [38[ the chirally improved Dirac 
operator 84| was used to evaluate this in the simplified local form 

F; v {x) ex ]TA 2 ^0r,*|A), 



A 

i 
2 



with the color generator T a in the fundamental representation. Zero modes do not contribute 
to this representation of the field strength. The ultraviolet filtering consists now in including 
only a certain number of low-lying modes in this spectral representation. 

Here we explore the application of this formula to the eigenmodes of the overlap Dirac 
operator -D(O) and compare the revealed structure with that shown by the ultraviolet filtered 
topological density. Obviously, a normalization factor remains undetermined, in particular 
because only a subset of low-lying modes will be included in the filter. 



In two recent papers 39|, |40|], a representation of field strength and action density has 
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FIG. 23: Normalized histograms with respect to the local (anti-)selfduality of the field strength 
tensor in the Q = subsample consisting of 37 configurations generated on the 16 3 x 32 lattice at 
(3 = 8.45, (a) taken over 100 % of the lattice sites, but depending on the number of non-zero modes 
(2 - 20) included in the filter, (b) applying various cuts specifying the infrared action density for 
the case of 20 eigenmodes included in the filter. 

been discussed, in a way analogous to Eqs. (Tl8l) and (fl9l) . which give the representation of 
the topological charge density directly in terms of the overlap operator. 



B. Local (anti-)selfduality of the infrared field strength tensor 

We now study the contribution of individual non-zero modes to the filtered field strength 
tensor. The normalization of the infrared field strength F" (x) does not come automatically. 
It could be obtained requiring that the infrared topological charge density 

qm(x) oc Tr (^(a;) F^x)) (36) 

with the dual field strength = ^e^upa F pa fits the ultraviolet filtered topological density 
q(x)x cut (for a suitable cut-off A cut ). The latter can be directly derived from the overlap 
operator following Eq. (fl9|) . A less elaborate normalization for Q ^ configurations could 
be that the total topological charge obtained from Tr (^F^ v (x) F fiu (x)j should reproduce the 
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index Q of -D(O) for the given configuration. 



One can define analogously the infrared action density 



sjr{x) oc Tr(F^ u (x) F^x)) 



(37) 



in order to compare it with q~m(x). With both infrared densities at hand, one can ignore the 
quest for normalization. It is now possible to turn the interest to the local (anti-)selfduality 
based on a ratio similar to Eq. (|30|) . 



analogously to Eq. (I3T1) . 

In Fig. [23] histograms are presented with respect to R(x) applying various cuts specifying 
the local action density §ir(o;). Firstly, sharp peaks at R = ±1 are observed that become 
weaker with the inclusion of more low-lying pairs of non-zero modes. We notice that even 
with 10 pairs included these peaks are still visible. They become more pronounced again 
when one focusses on part of the lattice sites applying a cut with respect to the infrared 
action density 5ir(x). 

C. Cluster analysis of selfdual and antiselfdual domains 

Since the close neighborhood of R = ±1 of these histograms corresponds to lattice points 
where the infrared field strength tensor is nearly ( ant i-) selfdual, it is interesting to see 
whether they are completely disconnected or form connected regions in space-time. We 
construct isosurface plots of |-R(x)| = -R cut = 1 — e demonstrated in Fig. [24] for a filter 
including the lowest 5 pairs of non-zero modes. The bubbles denote regions where the 
infrared field strength tensor is more (anti-) selfdual inside than outside. Red refers to 
antiselfdual regions, green to selfdual ones. Similar clusters have been seen when we have 
plotted clusters of the ultraviolet filtered topological charge density. 

In order to get a more quantitative characterization of the vacuum, it is useful to perform 
a cluster analysis on the basis of the quantity R(x). In Fig. [25] we show how the cluster 



sm{x) - q m {x 
sm{x) + q m {x 



(38) 



This is converted to R(x) 




(39) 
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FIG. 24: Isosurface plots of R(x) in the order of decreasing tolerance with respect to deviations from 
perfect (anti-)selfduality, (a) R cut = 0.97, (b) R cut = 0.98, (c) R cut = 0.99 and (d) R cut = 0.999. 
The figure shows the timeslice t = 6 of the typical Q = configuration on the 16 3 x 32 lattice 
generated at (3 = 8.45. Red and green surfaces enclose regions where antiselfduality or selfduality 
is better fulfilled than in the rest of the volume. The filter is using the 5 lowest pairs of non-zero 
eigenmodes. 
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FIG. 25: i?-cluster analysis using a running cut-off R C ut on the basis of the filtered action density 
sir and topological charge density used to define R(x). The i? C ut dependence is shown of (a) 
the total number of separate clusters, (b) the size of the largest cluster relative to all clusters, 
(c) the distance between the two largest clusters and (d) the connectivity (see text). The data is 
averaged over the 16 3 x 32 lattice ensemble generated at /3 = 8.45. 

composition changes with varying -R cut . When percolation begins close to R cut = 0.98 the 
number of separate .R-clusters starts to slightly decrease with decreasing -R cut . This means 
that allowing more tolerance with respect to violation of ( ant i-) self duality inside the clusters 
leads to the growth and coalescence of already identified clusters, rather than to the detection 
of further, isolated clusters. One can see that the distance between the two largest clusters 
remains relatively large (decreasing only from 12 to 7), meaning that a lot of space-time 
remains unoccupied by the almost (anti-)selfdual -R-clusters, even in the percolating regime 
signaled by the connectivity becoming f(r max ) ^ 0. The unoccupied regions are regions 
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FIG. 26: Results of the .R-cluster analysis at R cut = 0.96 (left) and i? cut = 0.99 (right) concerning 
the relative volume of all emerging clusters ((a) and (b)) and concerning the charge inside the 
clusters ((c) and (d)). The cluster charge is defined by the all-scale topological charge density 
(|18p . The distributions describe the same Q = configuration from the 16 3 x 32 lattice ensemble 
generated at j3 = 8.45 already considered in Fig. [6l and [241 

where the infrared field strength tensor is violating selfduality or antiselfduality more than 
the given tolerance -R cut permits. 

It turns out again that one largest cluster exists, which is rather large already before 
the onset of percolation, occupying 40 % of the total volume of all clusters. This i?-cluster 
rapidly becomes the absolutely dominating one with respect to the volume (more than 50 
%) after percolation has set in. This fact (topological "clumping") remains unexplained 
in all instanton/caloron models which are based on uncorrelated (anti-)selfdual objects of 
Q — ±1, eventually taking the attractive correlation of opposite charge cluster into account 
(that is actually seen at higher resolution A cut ). 

How the whole cluster collection changes between the percolating (R cn t = 0.96) and the 
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(c) (d) 

FIG. 27: Direct comparison between the isosurface plots with respect to R (left: with i? cu t = 0.99 
in (a) and i? cut = 0.98 in (c)) and the isosurface plots with respect to the truncated topological 
density q\ cut with A cut = 200 MeV (9 pairs of non-zero modes) (right, with <7cut/<2max = 0.2 in (b) 
and (/cut/<?max = 0.1 in (d)). All subfigures show the same time slice t = 6 of the typical Q = 
configuration of the 16 3 x 32 ensemble at /3 = 8.45. The similarity is found in all timeslices. 
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(a) (b) (c) 




FIG. 28: Two-dimensional section through the typical Q = configuration generated at /3 = 8.45 
on a 16 3 x 32 lattice. Direct comparison of (a) qm(x) based on the infrared field strength tensor 
based on 10 non-zero modes (5 pairs), (b) the mode-truncated density q\ cut {x) with A cut = 200 MeV 
(for this configuration it amounts to 9 pairs of non-zero modes), (c) the all-scale charge density 
q(x). Subpanels (d) and (e) show the positive and negative part of the local (anti-) self duality 
variable R(x). Note the 20-fold larger vertical scale in (c) compared to (b). The scale of (a) is 
freely adapted to that of (b). 

non-percolating regime (i? C ut = 0.99), with respect to charge volume, is illustrated in Fig. [26] 
for the Q = configuration already considered in Figs. [6] and [2H In order to define the 
cluster charge, the all-scale topological charge density q(x) f|T8|) has been summed over the 
sites belonging to the respective connected .R-clusters. 

Before percolation sets in (at -R cut = 0.99) one finds in this particular configuration O(20) 
significant clusters. Among them one cluster has an exceedingly large volume and two have 
a cluster charge Qciuster = 1 or —3. Deep in the percolation regime (at R cut = 0.96), there 
are essentially only two significant -R-clusters left over which share the total occupied volume 
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in the ratio 55/45. Their cluster charges are large and must compensate each other. Even 
in this state, the distance between the two clusters is still large, d(c\, C2) ~ 7 a. 

Near the maximum multiplicity, before the percolation sets in, for a suitable ultraviolet 
cut-off A cut the number of g-clusters is of the same order as the number of (anti-) self dual 
i?-clusters (domains) defined with a corresponding number of modes. Moreover, it turns 
out, that the respective isosurface plots are closely correlated. We consider this coincidence 
as an argument for the importance of locally (anti-) self dual infrared gauge fields as pert of 
the vacuum structure at a cut-off scale of 200 MeV. The correlation of positions and sizes 
of both types of clusters is shown in Fig. [27| where one timeslice is visualized by means of 
R(x) on the left side and of q\ cut (x) with A cut = 200 MeV (which actually corresponds to 9 
pairs of non-zero modes for this configuration) on the right side. The same kind of similarity 
exists in all timeslices. 

Summarizing, in Fig. [28] we show the various ways to present the topological structure in 
a two-dimensional section. The upper panels present the profiles of various densities, ?ir(x) 
based on the infrared field strength, the mode-truncated density q\ cut (x) and the all-scale 
density q(x). The lower panels show how the positive and negative part of R(x) restricted to 
the interval 0.95 < |-R(x)| < 1 highlights the regions of selfduality and antiselfduality. The 
plots illustrate this for the same Q = example configuration as visualized many times. 

VI. DISCUSSION AND CONCLUSIONS 

We have reported on several attempts to elucidate the structure of the (quenched) QCD 
vacuum using overlap fermions as a probe. In particular, we have looked at the localization 
properties of the eigenmodes and of the topological charge density which is built up by all 
modes. Furthermore, we have studied the chiral properties of the eigenmodes of the overlap 
Dirac operator and tried to relate this to the topological charge density. Finally we have 
searched for (anti-) self dual domains of the ultraviolet filtered field strength tensor. 

As a side-result of this work, we have obtained the topological susceptibility for the 
ensembles in use (from the index of each configuration) and the average spectral density. The 
analysis of the latter has resulted in an effective chiral condensate E e g for each lattice size. 
The final analysis of this in the light of quenched chiral perturbation theory, the comparison 
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with random matrix theory, in particular under the aspect of spectral correlations, will be 
subject of a separate publication. The emphasis of the present paper was put on the chiral 
and topological structural properties. 

While average properties as the dimensionality of eigenmodes (for zero modes and non- 
zero modes in special intervals of eigenvalues A) could be obtained by the scaling properties 
of the average inverse participation ratio (and its generalization based on higher moments 
of the scalar density) in agreement with critical spectral correlations (not described in this 
paper), the space-time structure of single modes and of the topological charge density could 
be studied by employing a cluster algorithm and, subsequently, a random walk algorithm. 
While the cluster algorithm locates and separates the clusters, depending on a cut (g cut or 
Pcut) applied to the quantity that is being considered, the random walkers help to specify 
the space-time dimension of the cluster by means of their return probability to the cluster 
center from which they start. 

A necessary condition for a pion to be able to propagate through a given lattice vacuum 
is that the scalar density of the modes forms percolating clusters, such that any two separate 
points on the lattice must be "connected" (in the sense defined in Section Hi El with non-zero 
measure when integrated over all configurations). Similar properties might be established 
for other quantities, like the topological charge density. 

In the specific case of the topological charge density one finds at low enough density 
only two connected clusters of opposite charge, which are infinitely thin in the co-direction 
and form a multi-layered compound that covers all space-time. This result can hardly be 
reconciled with the picture of a dilute gas of semi-classical, finite-sized, in four dimensions 
coherently extended excitations (instantons or other classical solutions). Since the topolog- 
ical charge density correlator in the vacuum must be negative for all non-zero distances (on 
the lattice at least beyond a distance of 0(a) set by the non-ultralocal density definition in 
terms of the overlap operator), this requires a local fractal dimension d* < 4 for the topo- 
logical density. The simultaneous presence of infrared (long-range) and ultraviolet (0(a)) 
features motivates to call the overlap topological density "all-scale density", as long as no 
mode-truncation is applied. This picture is known from the work of Horvath et ai, but it 
does not exclude the presence of even lower dimensional structure inside the two clusters. 
Here we complete the picture by pointing out (see Table V in the Appendix) that percolating 
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FIG. 29: The full pion propagator compared with the contribution of the zero modes alone and 
the cumulative approximation by various numbers of non-zero modes (see legend box) is shown in 
analogy to Fig. [TJ For this calculation 250 configurations have been used from the 16 3 x 32 lattice 
ensemble generated at j3 = 8.45. 

clusters of topological charge already appear with lower dimensions d* — 1 — 2 at higher 
q cut « 0.25g 

max • 

Similar things are observed for single modes of the overlap operator. Percolation sets in 
at dimensions d* ~ 2 and p cnt ~ 0.20p max for zero modes and the first few non-zero modes, 
whereas higher modes percolate already at higher p C ut/Pmax- 

In the Introduction (in Fig. [I]) we have seen that the long-distance part of the pion 
propagator is well saturated by the lowest 40 eigenmodes of the overlap operator. In more 
detail, Fig. [291 shows the pion correlator separately for an increasing number of modes, from 
the zero modes alone up to the additional inclusion of 40 non-zero modes. 

In Fig. [30] we show the convergence to the full propagator separately for various time 
differences t and in Fig. ED the effective pion mass separately for the full propagator and 
different mode-truncated approximations to it, also as a function of t. We see that the pion 
propagates with approximately 70 % of its physical mass even when the quark propagator 
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FIG. 30: The dependence of the pion propagator for various t (see legend box) on the number of 
accumulated non-zero modes is presented for the 16 3 x 32 lattice ensemble generated at = 8.45. 

is stripped down to the zero modes. 

If we wanted to define the mobility edge by the lowest mode at which propagation sets in, 
we would have to conclude that the mobility edge is at A = 0, and that the (quenched) QCD 
vacuum is critical in the whole lower part (A < 200MeV) of the spectrum, very analogous to 



871 ] . The level compressibility (not 



certain 2D or 3D disordered condensed matter systems 
described in this paper) supports this analogy. 

What changes in the vacuum picture if we would restrict our analysis of the vacuum 
structure to the lowest part of the spectrum? If the topological density is ultraviolet filtered 
by truncating its definition to the lowest modes, the emerging picture is completely different 
from what we have outlined above. The topological charge density correlator would now 
be resembling the presence of clusters of finite extent being subject to some weak but non- 
trivial intercluster correlations, increasing with the cut-off A cut . This picture is basically 
reproduced by the cluster analysis and the random walker algorithm. 

In the Appendix we summarize our main results concerning the dimensionality of the 
all-scale and of the mode-truncated (ultraviolet filtered) topological charge density (Tables 
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FIG. 31: The effective-mass plot vs. t for the 16 3 x 32 lattice ensemble generated at j3 = 8.45. 
This is shown for the full pion propagator and for various numbers (see legend box) of non-zero 
modes included in the approximation. 

V and VI) and of the scalar density of zero modes and selected non-zero modes (Tables 
VII and VIII). Here we also test the reliability of the "covering sphere method" introduced 
in subsection IIVDI that attempts to estimate the effective dimensionality by the growth 
of the integrated quantity under consideration (the topological charge and scalar density) 
with the radius of a sphere covering a cluster. Although this latter method clearly defines 
a certain cut-off that separates a linear growth from a nonlinear one in both cases, the 
estimated dimensionalities do not coincide with the random walker method. This shows 
that the density inside the clusters is very inhomogeneous. 

Concerning the all-scale topological density, it is remarkable that at low values of the 
density (at the lowest cut-off q cu t defining the clusters) the dimensionality is practically 
independent of the lattice cut-off a and becomes almost equal three. The covering sphere 
method finds it somewhat bigger than three. At higher values of the density smaller dimen- 
sionalities are found. There the covering sphere method gives an effective dimension equal 
to one, i.e. the charge inside a cluster grows linearly with the size. For the mode-truncated 
density at lowest values of the density (at a low cut-off defining the clusters) the clusters' 
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dimensionality surprisingly is also close to three. How the infrared clusters (with a fixed 
A cut ) start to percolate we can describe in the following way: at high values of the density 
they are well-separated and localized to zero dimensional balls or one dimensional rods of 
finite thickness, before they start to percolate at somewhat lower density with an effective 
dimension still relatively low, d* = 1.5 — 2. Finally, at even lower values of the density, 
they become extended in essentially three dimensions with an a-independent thickness of 
O(Aqcd). 

A nonvanishing local chirality X\(x) is assumed by individual lower-lying non-zero modes 
in regions of their preferential localization (i.e. at high scalar density p(x)). This effect 
could be demonstrated to be correlated to the ultraviolet filtered topological density q\ over 
a distance comparable to the correlation length of the density itself. Here we can only say 
that this ultraviolet filtered density is resulting from the correlated behavior of all lowest- 
lying modes. On the other hand, this type of correlation of X\ with the topological charge 
density is a natural consequence in the instanton model where the topological charge is 
simply carried by the semiclassical background and the non-zero modes are superpositions 
of the (anti-)instantons' zero modes. 

The all-scale topological density does not follow this pattern. It deserves its name at 
low cut-off g cut where it combines long-range coherence on one hand and lower-dimensional, 
i. e. singular structure on the other. When the cut-off is increased, percolation breaks down 
and an increasing number of isolated clusters appears with lower and lower dimensions. In 
this sense the all-scale topological density is multifractal, too. It should be noted that the 
maximal number of these clusters seems to depend only on the lattice size in lattice units. 
In a given physical volume it is diverging with the lattice spacing a — > 0. 

Applying the same methods to the scalar density of individual modes (from zero modes 
up to the highest analyzed modes) we also have found that at very low cut-off p cut they do 
not extend isotropically but show a lower-dimensional structure. At the lowest cut-off the 
effective dimension ranges from two-dimensional for zero modes and the very first non-zero 
modes to three-dimensional for the modes up to the 120th non-zero mode. This means that 
all the modes that are sufficient to saturate the pion propagator share this property of low- 
dimensionality. This is in agreement with the critical level statistics (being "non-metallic" 
following the condensed-matter terminology). 
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Summarizing we can state that at sufficiently low density (q(x) or p(x)) we find consis- 
tently a dimensionality d* ~ 3 to characterize most considered densities. The only exception 
are the very lowest modes which show a two-dimensional structure with respect to the growth 
of the integrated norm. This special feature is corroborating the average properties that have 
been concluded from the scaling behavior of the IPR. 

It is remarkable that the mode truncated topological charge density percolates already at 
relatively low (volume) packing fraction, when its domains have a finite extension, forming 
coherent clusters of finite thickness O(Acd)- This is impossible for the all-scale topological 
density. This kind of clustering of the ultraviolet filtered topological charge is corroborated 
by the search for ( ant i-) self dual domains, employing the ultraviolet filtered field strength 
tensor. The latter can be defined by means of the overlap operator, too, employing filtering 
in the form of truncating the spectral representation. 

In conclusion: the vacuum picture exhibited by the modes that are sufficient to build 
up the pion propagator resembles a dilute liquid of separated ( ant i-) self dual excitations if 
looked at it at higher cut-off g cut of the topological charge density and a dilute, percolating 
network at lower cut-off. The less-than-three dimensional pinning-down of the very lowest 
modes at very high scalar density p(x) (and of the topological density at very high values of 
q\{x)) might be related to singular (confining) excitations like monopoles or vortices. This 
possibility is presently under investigation. In particular, the two-dimensional structure of 
the zero modes (with pointlike peaks on top of the latter) intriguingly suggests that they 



are pinned to vortex surfaces and vortex intersections [85(. In this sense we may conjecture 
that the lowest part of the fermion spectrum is both realizing chiral symmetry breaking and 
tracing the agents of confinement. This would be close in spirit to recent attempts in the 



literature to extract t 



Dirac operator 86 



88 



re dynamics of the Polyakov loop from the spectral properties of the 

m. 



It will be our next task to extend these studies to finite temperature and full QCD in 
order to disclose the changes of topological structure under the influence of dynamical quarks 
and those accompanying the chiral phase transition. 
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VII. APPENDIX: TABLES OF DIMENSIONS 



Qcut/lmax 


0.00 0.10 0.20 0.30 0.40 0.50 


12 3 x 24 all-scale 
12 3 x 24 A cut < 200 MeV 
12 3 x 24 A cut < 800 MeV 


2.7(1) 1.5(1) 1.1(1) 0.7(1) 0.4(1) 0.1(1) 
3.1(1) 1.8(1) 1.5(1) 1.0(1) 0.7(1) 0.2(1) 
3.0(1) 1.8(1) 1.5(1) 1.0(1) 0.7(1) 0.2(1) 


16 3 x 32 all-scale 
16 3 x 32 A cut < 200 MeV 
16 3 x 32 A cut < 634 MeV 


2.6(2) 1.5(1) 1.1(1) 0.9(1) 0.5(1) 0.1(1) 
3.0(1) 2.1(2) 1.8(1) 1.2(1) 0.8(1) 0.3(1) 
3.0(1) 1.9(1) 1.7(1) 1.2(1) 0.7(1) 0.3(1) 



TABLE V: Effective dimension d* of the leading cluster of the topological charge density based on 
the random walker approach. Dimensions referring to the percolating regime are printed in bold 
face. The data is averaged over the 12 3 x 24 (16 3 x 32) lattice ensemble at (3 = 8.10 (f3 = 8.45). 



9cut/?max 


0.00 0.10 0.20 0.30 0.40 0.50 


12 3 x 24 all-scale 
12 3 x 24 A cut < 200 MeV 
12 3 x 24 A cut < 800 MeV 


3.3(1) 3.1(1) 1111 
2.5(2) 11111 
2.9(1) 1.8(1) 1111 


16 3 x 32 all-scale 
16 3 x 32 A cut < 200 MeV 
16 3 x 32 A cut < 634 MeV 


3.3(1) 3.3(1) 2.5(1) 1 1 1 
1.9(1) 11111 
2.3(1) 11111 



TABLE VI: Effective dimension of the leading cluster of the topological charge density based on 
the covering sphere method. The data is averaged over the 12 3 x 24 (16 3 x 32) lattice ensemble at 
/? = 8.10 (0 = 8.45). 
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Pcut/Pmax 


0.00 


0.10 


0.20 


0.30 


0.40 


zeromodes 


3.6(2) 


2.8(2) 


2.0(2) 


1.4(1) 


1.1(1) 


1st non-zero mode 


3.5(1) 


2.9(2) 


1.9(2) 


1.4(1) 


1.1(1) 


10th non-zero mode 


3.6(2) 


3.0(2) 


2.2(2) 


1.6(1) 


1.1(1) 


30th non-zero mode 


3.5(1) 


3.2(2) 


2.6(1) 


2.0(1) 


1.4(1) 


50th non-zero mode 


3.6(2) 


3.3(1) 


2.8(1) 


2.1(1) 


1.5(1) 


70th non-zero mode 


3.6(2) 


3.4(1) 


2.9(1) 


2.2(1) 


1.6(1) 


90th non-zero mode 


3.5(2) 


3.3(2) 


3.0(1) 


2.4(1) 


1.7(1) 


120th non-zero mode 


3.6(2) 


3.4(2) 


3.1(1) 


2.5(1) 


1.8(1) 



TABLE VII: Effective dimension d* of the leading cluster of the scalar density based on the random 
walker approach. The data is averaged over the 16 3 x 32 lattice ensemble generated at (3 = 8.45. 
Dimensions referring to the percolating regime are printed in bold face. 



Pcut/Pmax 


0.00 


0.10 


0.20 


0.30 


0.40 


zeromodes 


1.9(1) 


1 


1 


1 


1 


1st non-zero mode 


2.1(1) 


1 


1 


1 


1 


10th non-zero mode 


2.5(1) 


1 


1 


1 


1 


30th non-zero mode 


2.9(1) 


2.3(1) 


1 


1 


1 


50th non-zero mode 


2.9(1) 


2.4(1) 


1 


1 


1 


70th non-zero mode 


3.0(1) 


2.7(1) 


1.9(1) 


1 


1 


90th non-zero mode 


3.1(1) 


2.6(1) 


1.9(1) 


1 


1 


120th non-zero mode 


3.1(1) 


2.9(1) 


2.3(1) 


1 


1 



TABLE VIII: Effective dimension of the leading cluster of the scalar density based on the covering 
sphere method. The data is averaged over the 16 3 x 32 lattice ensemble generated at (3 = 8.45. 
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